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Abstract 

A formalism for studying heavy quarks in terms of model Dyson-Schwinger 
equations is developed. The formalism is the natural extension of a technique 
which has proved successful in a number of studies of light hadron physics. 
The dressed heavy quark propagator, calculated to leading order in the in- 
verse quark mass, is incorporated in a treatment of mesons consisting of a 
heavy quark and light antiquark via the ladder approximation Bethe-Salpeter 
equation. In the limit of infinite heavy quark mass the model is found to 
respect the spectrum degeneracies present in Heavy Quark Effective Theory. 
An exploratory numerical analysis of a simple form of the model is carried 
out to assess its viability for studying D and B mesons. 
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I. INTRODUCTION 



There has been a great deal of recent theoretical interest in heavy quark mesons, that is, 
mesons containing charm quarks or bottom quarks such as B and D mesons. The basis of 
this interest is heavy quark effective theory (HQET) HJ, which hinges on the realisation that 
in the limit of infinite quark mass extra symmetries, and hence degeneracies, occur beyond 
those normally associated with QCD 0. It is a property of many existing HQET analyses, 
particularly of meson spectrum calculations, that dynamical self dressing of the quarks is ei- 
ther not included at all, or rarely included in a non-perturbative fashion. Examples of recent 
spectrum calculations in which heavy quarks are modelled essentially by free propagators 
are given in ref . @ . However, one could argue that a theory which does not prevent quarks 
from being on mass shell is not in keeping with the spirit of a confinement and may lead, in 
certain instances, to incorrect results. In Bethe-Salpeter equation calculations for instance, 
there is in principle nothing to prevent the occurrence of spurious production thresholds. 

At the other end of the meson spectrum, namely the physics of light quarks, the Dyson- 
Schwinger equation (DSE) technique j| has enjoyed considerable success in modelling a 
number of light meson phenomena. Basic to the DSE technique is the quark propagator. 
The spacelike structure of this propagator has been well modelled within the DSE technique 
for light quarks. At a phenomenological level, model quark propagators based on numerical 
studies of approximate Dyson Schwinger equations have proved efficacious in calculating light 
meson spectra 0-0, electromagnetic form factors of the pion and kaon |5J, tt-ti scattering 



lengths ||, quark loop contributions to p-u mixing |10[ and the anomalous •yn — > 7 and 
77r — > 7T7T form factors |TT|| . 

The aim of this paper is to formulate an extension of the DSE technique to the realm 
of heavy quarks. In the applications of the the DSE technique to light quarks listed above, 
it has not been necessary to determine the behaviour of quark propagators well away from 
the real spacelike p 2 -axis. However, for the purpose of accurately modelling confined heavy 
quarks, because the self energy is small compared with the current quark mass rriQ, it is 
particularly important to know the analytic structure of the heavy quark propagator near the 
(timelike) bare fermion mass pole ^Minkowski = m Q- The method frequently used within the 
DSE technique of fitting analytic functions to the numerically determined spacelike quark 
propagator is inadequate for this task. It is for this reason that new ideas must be developed 
within the DSE technique if it is to be adapted to heavy quarks. 

It is the philosophy of the DSE technique that confinement in QCD should be signalled 
by the absence of timelike poles in the quark propagator. This analytic structure is the 
result of non-perturbative quark self dressing, which in principle could be demonstrated by 
a full treatment of the quark Dyson-Schwinger equation. Many numerical calculations of 
approximate quark Dyson-Schwinger equations in confining field theories have demonstrated 
propagator solutions which are free of timelike poles M, though not necessarily free of 



complex conjugate poles away from the real p 2 -axis |TI|. Herein we shall examine in detail 



the analytic structure of the heavy quark propagator in the limit of infinite quark mass. We 
shall also consider the Bethe-Salpeter equation (BSE) for bound states of heavy quarks and 
light antiquarks. Our eventual aim is to make contact with the observed spectrum of D and 
B mesons. However, at the present level of approximation we are unable to find solutions 
to our model BSE because of spurious poles in the model quark propagators. We suggest 
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ways in which the model can be improved to overcome this problem. 

The layout of the paper is as follows: In Section II we develop the basic formalism for the 
heavy quark propagator. In Section III we consider the Bethe-Salpeter equation for mesons 
consisting of one heavy and one light antiquark. Sections IV and V deal with scalar and 
vector mesons respectively, and demonstrate that the symmetries of HQET are preserved 
within the heavy quark limit of the DSE technique. Some numerical results are given in 
Section VI and conclusions drawn in Section VII. 

II. THE HEAVY QUARK PROPAGATOR. 

Before understanding the procedures which must be employed to ensure the extension of 
the DSE technique to heavy quarks, it is important to appreciate that a naive translation of 
the light quark sector DSE technique, by setting the quark mass to large values in the existing 
formalism and computer codes, simply will not work. This is because the interesting physical 
phenomena are driven by non-perturbative contributions to the heavy quark propagator in 
the vicinity of the bare fermion mass pole. These contributions are of order I/tuq compared 
with the bare fermion propagator. Earlier studies of the heavy fermion limit of another 
confining theory, QED3 [[OJ, demonstrate categorically that for fermions whose mass is 
greater than the typical scale of the interaction of the theory, this information cannot be 
accurately obtained by an extrapolation of the spacelike solution of the fermion propagator 
produced by the usual DSE technique, and a new formalism must be developed. 

We begin with the general form for the inverse of the dressed quark propagator consistent 
with Lorentz and CPT covariance: 

S-\p) =ii-vA(v 2 ) + B{ V 2 ), (2.1) 

where A and B are scalar functions. We use a Euclidean metric in which timelike 4-vectors 
satisfy p 2 < 0. Within the DSE technique the full quark Dyson-Schwinger equation is 
typically truncated to a manageable approximation. For the purposes of illustrating the 
extension to the heavy quark sector, we consider the simplest viable approximation to the 
quark DSE, 

4 f d 4 q 

S~ (p) -ij-p-m Q = - J — -—^(g^Afr - q), (2.2) 

in which the quark-gluon vertex has been replaced by a bare vertex, and the dressed gluon 
propagator is modelled by a representation g 2 D fJiU (k) = 5^ u ^.{k 2 ) in a Feynman-like gauge. 
Our formalism for calculating model heavy quark propagators can be extended without 
trouble to include more realistic vertex ansatze (to allow for gauge covariance |Q) and to 
more realistic model gluon propagators (see, for instance, @]). 

In the limit of heavy quark masses, tuq — > oo, the dressed quark propagator is dominated 
by the bare form S^ rc = ■ p + niQ. However it is important to isolate from the full 
propagator order 1 / mq corrections to the bare propagator which drive confining and remnant 
chiral symmetry breaking effects. To this end we set 

m = l + B ^ )=m J 1 + ^m, (2.3, 

m Q V m Q I 
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where, in the spirit of HQET, we have introduced a new momentum variable 



= irrigv^ + hp, (2.4) 

with t> M a constant unit 4-vector. 

The functions Ha arid X# can be expanded as a series in 1/itlq, and the DSE solved 
order by order as a set of coupled integral equations. In the appendix we derive the set 
of four coupled equations required to solve the heavy quark propagator to O^I/ttiq). To 
leading order, we find that and E# can be written as functions of the single variable k-v, 
and the quark propagator becomes 

s & = 1+ o* -u iv/n v7i s + K 1 ) • ^ 

2 ik • v + Eb(« • v) — T,a,(k ■ v) v v/ 

Substituting into Eq. (|2.2j), changing the variable of integration to k' defined by q^ = imv fl + 
k', and projecting out coefficients of ff, ft and I gives a set of integral equations which can 
be summarised in the equation (see Eqs. ( |A.8| ) and ( |A.9[ )) 



where 



4 r d 4 k' 1 

E(k -v) — - / y—-rA(k — k')- =— -, (2.6) 

v J 3J 2tt 4 v J ik' ■ v + Z(k' ■ v) V ; 



E(k-v) =Z B {k-v) -E A (k-v). (2.7) 



Note that to leading order in 1/mq the heavy quark self energy is specified by a single, 
complex valued, scalar function. This equation is numerically tractable once the function 
A is specified, and provides the required information regarding fermion self dressing in the 
vicinity of the bare fermion mass pole, m 0, in the limit niq — > oo. 

Note also that the change of integration variable from q to k! involves an assumption that 
the propagator, and hence the functions A(q 2 )/[q 2 A(q 2 ) 2 + B(q 2 ) 2 } and B(q 2 )/[q 2 A(q 2 ) 2 + 
B(q 2 ) 2 }, be analytic over the region 



, , (Im(g 2 ))^ 



-m A Q + ± -^f ■ (2.8) 

In the limit iriq — > oo this translates to the requirement that l/[ik ■ v + E(/c • v)} be analytic 
for Im(k ■ v) < 0. 



III. THE BETHE-SALPETER EQUATION FOR HEAVY-LIGHT MESONS. 

The combination of bare vertex DSE and ladder approximation Bethe-Salpeter equation 
produces well the spectrum of light quark mesons Here we extend the formalism to 

mesons containing one heavy quark and one light antiquark in the limit mq — > oo. We 



note that our work differs in a fundamental way from that of Villate et al. JL5] who use an 
instantaneous approximation to study heavy quarkonium states via the combination of DSE 
and BSE. Because the light quark present in our case must be treated in a fully relativistic 
fashion, the instantaneous approximation is not applicable. 
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The homogeneous ladder approximation to the meson Bethe-Salpeter equation in a 
Feynman-like gauge is, 



r(p, P) = ~ J A(p - qftMq ~ ZP)T(q, P)S Q (q + (1 - t)P)^. (3-1) 

Here S q (p) is the full light quark propagator, Sq(p) the full heavy quark propagator, and 
S^Aip — q) the same effective gluon propagator as that employed in the DSE. The Bethe- 
Salpeter amplitude T(q, P) is defined so that external outgoing and incoming quark lines 
carry momenta q + (1 — £)P and q — £P respectively 
For the light quark, we write the propagator as 

s q(p) = ~i fav(p 2 ) + °s(p 2 ) 

= i p<A q (p 2 ) + B q (p*y (3,2) 

The vector and scalar parts oy and as are related to the functions A q and B q simply by 
dividing by the quantity p 2 A 2 (p 2 ) + B 2 (p 2 ). Using Eq. (|2.2| ), with the quark mass replaced 
by the light current quark mass m q , A q and B q are given by the integral equations 



p 2 (Mp 2 



d 4 q A q {q 2 ) 



16 /• d A q B q {q 2 



B ^ P ] ~ m « = t J w^ A(p " g) - (3 - 4) 

For the heavy quark we use the leading order heavy quark propagator Eq. (p3|). 



S Q (im Q v, + kfj) = l A^ lk . v + nk . v) (S-5) 

where E(fe • u) is obtained from Eq. ( p.6|) . 

We choose to work in the rest frame of the meson and set 

P„ = (0,i(m Q + 5)) (3.6) 

where 5 is the contribution to the meson mass from the "brown muck" of light quarks and 
glue surrounding the heavy quark. Our aim is to obtain the BSE in the limit tuq — ► oo in a 
form containing 5, but independent of tuq. It is convenient to reparameterise the momentum 
partitioning in Eq. (|3.1|) via 

(=^TJ. (3-7) 
rriQ + 5 



Since the light quark Eqs. (3.3) and ( |3.4|) are solved in the first instance for real, spacelike 
momenta, we see that the choice 77 = has the advantage that no numerical extrapolation 
of S q to complex momenta is required. On the other hand, using the freedom inherent in 
to set Vfj, = (0, 1), we see the choice r] = 1 has the advantage that no extrapolation of 
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the numerical solution of the heavy quark self energy £(^4) is required away from the real 
&4 axis. In practice, one expects that the optimum choice of 77 requiring minimal numerical 
extrapolation of either propagator lies somewhere between these values. 

Before deriving integral equations for scalar and vector heavy-light mesons, it is conve- 
nient make some definitions. Define the projection operators 

A± = (l± 74 )/2. (3.8) 

The light quark propagator can then be written as 

S q (q - £P) = S+A + + S_A_ + S ± ft (3.9) 

where 

H = (?iTi + 9272 + qsis)/ |q| , (3.10) 



with |q| = \Jq\ + q\ + q\ and 

^+(^4, |q|) = <rs(s) - {r}6 + iq 4 )a v (s), (3.11) 

S-(qt, |q|) = °s(s) + (r)6 + iq 4 )a v (s), (3.12) 

Sj_(<?4, |q|) = -i |q| crv(s), (3.13) 

where 

s = |q| 2 + q\ - r] 2 S 2 - 2ir]5q 4 . (3.14) 
With = (0, 1), the heavy quark propagator is, from Eq. ( |3.5| ), 

S Q (q + (1 - g)P) = , A + — . (3.15) 

m - {l-ri)6 + H (q 4 + i(l - 7])6) 

We see that, by judicious choice of momentum partitioning, the quark propagators, and 
hence the BSE, are independent of mq in the tuq — > 00 limit. 

IV. SCALAR AND PSEUDOSCALAR MESONS. 

The general forms of the pseudoscalar and scalar Bethe-Salpeter amplitudes are given 
by0 



r pseud (g, P) = (fi + hJP + h1+ Ml P}) 75, (4.i) 



T^(q, P) = f 1 + f 2 P + f^+ f A [l p], (4.2) 



where fi to / 4 are functions of q , P and q ■ P. With the P^ chosen as Eq. ( |3.6| ), these can 
equally well be written in terms of the projection operators A± as 
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r pseud (g) = (/+A+ + ; _ A _ + g+ + g _ 



(4.3) 



pcala r(g) = /+A+ + f A _ _ g+ ^ + _ g _ ^ (4 4) 

where f± and g± are functions of and |q|. 

By substituting Eqs. ( |3.9|) and ( |3.15|) for the quark propagators into Eq. ( |3.1| ) and making 



use of the identities 

7m A ±7m = A± + 3A T , 7m flA ±T/1 = -&. (4.5) 

a set of coupled integral equations can be projected out. For the pseudoscalar meson we 
obtain 

J - W 37 (2tt)^ IP q) iq A -(l- r,)S + E(q, + i(l — 77)^' 1 J 

/+(p) = 3/_(p) (4.7) 

/ \ 4 /■ A _ S+g-(q) + S ± f-(q) 

3 7 (27r) 4 zg 4 - (1 - ?7)d + S (54 + z(l - rj)5) 

9+(p)=9-(p)- (4-9) 

The scalar equations can be obtained from these by making the replacements as — > —<?s in 
Eqs. |U1|) and ( pTT2|) . 

V. VECTOR AND AXIAL VECTOR MESONS. 



Llewellyn-Smith |Tj| gives the general form of the J — 1 BS amplitude as 



where 



T,(q,P)=T,(A + q-PB 4 + C P + D[l P}) 
-Q^q-PB + 2 PD) 

+Q,(E + F 4+q-PG P + H[l P]) , (5.1) 



and A to H are functions of q 2 , P 2 and q ■ P. The second line of Eq. (|5.1| ) serves the purpose 
of ensuring that A to H are even functions of q ■ P if the fermion is odd under charge 
conjugation. Since we are interested in q-Q mesons, charge conjugation is not important, 
and this line can be ignored. 

Working in the vector meson rest frame, Eq. ( |3.6| ), the general form of the amplitude 
can be recast in the form = (T, 0) with 
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f(q) = f+(g)A+ + f-(g)A- 



(5.3) 



where 



(5.4) 



with /±, g±, /i± and /c± functions of 94 and |q|, and 7 = (71,72,73) and q = (91,92,93)/ |q|- 
Using the identity 7 q 7 = q + ~ [7, q] , it is convenient to absorb part of the g± terms into a 
redefined h± and write instead 



(5.5) 



The following relations are then useful for extracting the scalar coefficient functions from 



/±(p) = \ {Sij - M') tr (l^t 



(5.6) 



»±(p) = -igtr (h,m(p) 



(5.7) 



(5.8) 



fc±(p) = o (3p»Pj - 5y) tr (7 fc rf (p)) . 



(5.9) 



Substituting Eq. ( p~5|) into the BSE Eq. fl3.1|) together with the propagators Eqs. ( pT9| 
and ( |3.15| ), and making use of the identities 



7/^A±7m = A± + 3A T , 7^7iA±7 At = ~7i(A± + A T ) 



(5.10) 



7 M [ 7i , q]A ±7/ , = [ 7i , q](A± — A T ), 
we obtain the following set of coupled integral equations: 

4 f d 4 q 



U(p) 



-A(p-9)X 



3 J (2vr) z 

+ - ~ (S^k + + S ± (^ + + h+)) (l 

1 



(p-q) : 



iq 4 - (1 - ?7)5 + £ (94 77)5) 



/-(p) = /+(p), 



x 



(5.11) 



(5.12) 
(5.13) 
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9 + {p) = ~J (04 A (^-^)x 
(S + g + (q)-S ± f + (q))p.^ 



iq 4 — (1 — r})5 + E (c/4 + i(l — rj)8) ' 
0-(p) = -ff+(p), 



M ^i/(4 A(p_9)x 

(S + h + (q) + S ± (f + (q)+k + (q))pAr 



iq 4 -(l-r))6 + Z(q 4 + i(l-r))5y 
h_(p) = 3h+(p), 



\ [S^k + (q) + S L (g + (q) + h + (q))} (l - 3(p.q) 



iq 4 - (1 - rj)S + E (q 4 + i(l - 77)^) 



= fc+(p). 

These equations can be simplified greatly by setting 

f = 3f + + k + , 

g = h + - 2g + , 

~h = g+ + h+, 

h — . 

This gives 

f(v) = A [ d ' q A(v a) S-^M + S^m 

J[P) 3 J (2tt) 4 [P q >iq A -{\- r})S + E(q 4 + i(l-rj)Sy 

W 3 7 (27T) 4 P ' qZAlP ^-(l-^ + E^ + ^l-T^)' 

fc(p) = -it*3-i 6A( P - q) S + (q)h(q) + S ± (q)k(q) 

n[P) 3 J (2ny P - qlX[P q) tq4 -(l- v )5 + X(q 4 + t(l- V )5y 



Hp) 



d 4 q 



l-3(p.q) 2 )A(p-g) 



S_(q)k(g) + S ± (q)h(g) 
(l- V )5 + Z( qi + i(l-r))5y 



(5.27) 



Eqs. ( |5.24|) and ( |5.25| ) are equivalent to the scalar equations Eqs. fl4.6|) and (|4.8| ). They 
represent the 1~ meson whose light quark spin content is j = |, which, as predicted by 
HQET, is degenerate with the scalar 0~ meson. Eqs. ( |5.26j ) and ( |5.27|) represent the Y~ 
meson whose light quark spin content is j — |. 

The 1 + meson Bethe-Salpeter amplitude is obtained by multiplying the general form 
Eq. ( |5.3| ) through by 75. The corresponding set of integral equations can be obtained from 
those above by making the replacement as — > —<7s in Eqs. ( |3.11|) and ( |3.12|) . 



VI. NUMERICAL CALCULATIONS. 



In order to explore the potential of our formalism we consider here the simple gaussian 
gluon propagator ansatz 



A(p - q) 



3 (27r) 4_^ e -(p- 9 )V«. 

16 Oi 7T 



(6.1) 



This model propagator is simple enough to render the angular integrals in both the SDE and 
the BSE analytically tractable, but has the disadvantage that it underestimates the gluon 
interaction strength in the ultraviolet. However, the evidence from light meson applications 
of the SDE technique is that it is the infrared behaviour of the gluon propagator which most 
influences the low energy physics. Furthermore, the large mass expansion makes sense for 
non-perturbative interactions, but not for hard gluons. By using a gaussian model gluon 
propagator we provide an ultraviolet cutoff within a treatment which does not admit hard 
gluon contributions. We therefore regard this model propagator as an acceptable starting 
point. For small a Eq. (|6.1|) reduces to the infrared dominant model 



lim A(p — q) 



16 



(27r)V5 4 G9-g), 



(6.2) 



which has frequently been employed in exploratory applications of the DSE technique to 
light quarks 

Substituting Eq. ( p.l|) into the heavy quark DSE Eq. ( |2.6|) with v M = (0, 1) gives 



S(fc 4 



4(a7r) 1 /2 Jo 



elk a 



Here we have assumed £(— k^) = ^(k^)* 
negative k^ of the solution 



+ 



e -(k 4 +k' 4 f/a 



(6.3) 



ik' 4 + X(ki) -ik' A + Y,(k l A )\ 
this being a property of the obvious extension to 



lim E(&4) 



-ik 4 + J fi 



/04 



if < A; 4 < ji, 



(6.4) 



-k 4 + 



k~f 



if 



h > 



10 



to the algebraic equation arising from the infrared dominant propagator Eq. ( |6.2| ). Sub- 
stituting the gaussian propagator Eq. (|Q| ) into the light quark DSE Eqs. ( |3~3l) and ( |3.4|) 

gives 



Mp 2 



ap" Jo 



q 3 dq 



Mq 2 



-(p 2 +q [ 



q 2 A q (q^ + B q (q 



2\2 



(6.5) 



where I\ and I2 are modified Bessel functions. 

Here we shall carry out an initial study of the viability of the model gaussian gluon 
propagator in the light quark chiral limit m q — > 0, and heavy quark limit niq — > 00. In this 
limit the parameter a in Eq. ( |6.1| ) can be scaled to unity in all equations by making the 
replacements 

fx = y/otfi, 5 = y/aS, (6-7) 



E(h) = yfa£(h), A( v 2 )=A(f), B(p 2 ) = V^B(p 2 ), (6.8) 

where 

k^ — k^/y/a, p = p/\fa. (6.9) 

It is then sufficient to explore the parameter space of the single parameter fx. 

We have solved Eqs. fl6.5|) and ( |6.6| ) by numerical iteration for real, spacelike p 2 for a 
range of values of fx. To extrapolate the light quark propagator from the spacelike axis into 
the complex p 2 -plane, the right hand sides of Eqs. (|3~3l) and (|3.4j ) were integrated along the 
real spacelike g 2 -axis using the converged solutions together with complex values of p 2 . We 
find that for fx < 3, the light quark propagator has a timelike pole, and that for larger values 
of fx the pole moves away from the real axis into the complex plane. Table | lists the position 
p 2 = s polc of the pole as a function of fx. 

To obtain a numerical solution to Eq. ( |6.3j ) for real we iterated from Eq. ( |6.4| ). The 
numerical solution for fx = 3.5 is shown in Fig. 1. The function £ is extended to complex 
values of k 4 by integrating the right hand side of Eq. ( |6.3| ) along the real /c^-axis. We find 
numerically that the requirement stated at the end of Section III, viz. that the heavy quark 
propagator be free of singularities for Im/c 4 < 0, is satisfied over the range of fx listed in 
Table |. However, we also find a succession poles which migrate inwards from ioo along the 
positive imaginary fc4-axis before veering away from the axis as fx is increased. The position 
of these poles is plotted in Fig. 2, and the position of the pole closest to the real /c 4 -axis 
is listed in Table Q. We note also that, as fx —>■ 00 corresponding to the infrared dominant 
limit Eq. ( |6.2| ), the poles move further from the origin. This is consistent with the limiting 
solution Eq. flOl ) which is free from poles except at k± = 00. 

If the quarks are to be a confined particles, we are restricted to values of fx for which 
the light quark propagator S q (p) has no pole on the real timelike p 2 -axis, and for which the 
heavy quark propagator 1/ (ik^ + S(/c 4 )) has no poles on the imaginary fc 4 -axis. While this 
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requirement can be fulfilled by our light quark propagator solutions, unfortunately we find 
the heavy quark propagator is not totally free of timelike poles for any value of \x over the 
range considered. It is reasonable to suggest that this fault could be rectified by improving 
the approximations employed DSE (bare vertex and gaussian gluon propagtor). In fact, 
in ref. ||19|| , it is shown that if the infrared dominant gluon propagator Eq. ( |6.2| ) is used 
in conjunction with a vertex ansatz satisfying the Ward-Takahashi identity, one obtains a 
quark propagator which is an entire function of p 2 irrespective of the quark mass. It is 
therefore reasonable to assume that the timelike poles in the heavy quark propagator are an 
artifact of the model. 

Poles which occur away from the timelike axis are, in principle, no impediment to con- 
finement. However, solution of the BSE requires knowledge of the quark propagators over a 
region of the complex plane. From Eqs. ( |3.9| ) to ( |3.15| ), we see that the light quark propagator 
S g (p) is sampled by the BSE over the region 



Reg > (Img)2 -^ 

where s = p 2 , and the heavy quark propagator is sampled along the line 

Im fc 4 = (1 — r))S, 



(6.10) 



(6.11) 



where 5 is the "brown muck" contribution to the meson mass, and r] a momentum parti- 
tioning parameter. Assuming there are no compensating zeros in the BS amplitudes, the 
integral in the BSE is divergent if poles are encounted in either of these two regions. Writing 
-pole = x + iY and Im k\ oXe = K, it is straightforward to show that the model then fails if 



5 > 5„ 



K + 



VX 2 + Y 2 -X 



1/2 



(6.12) 



To approach the maximum allowed value of 5 in numerical calculations, it is necessary to 
choose the momentum partitioning rj to take its optimum value 



Vopt 



VX 2 + Y 2 -X 



1/2 



V2K + \^fX 2 TY 2 - X 



1/2- 



(6.13) 



Values of <5 max and rj opt are plotted in Table |. 

Solution of the BS equations involves the iteration of the sets of coupled linear integral 
equations obtained in Sections V and VI. For the J p = 0~ meson, for instance, from 
Eqs. (|4.6|) and Q4.8Q and the definitions ( |3.11|) to ( |3.13|) we can assume without loss of 
generality that 



/-(|q| , -il) = /-(|q| , 94)*, g-(h\ , -<&) = -#-(|q| , q^T- 



(6.14) 



Executing the spatial angular interations analytically then gives a set of equations of the 
form 



f(\p\ ,Pa) 



dq4 



d|q| K(\p\ ,p 4 ; |q| ,g 4 ;5)/(|q| ,g 4 ), 



(6.15) 
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where / = (Re/_, Im/L, Re g_, lmg_) T and the function K is a 4 x 4 matrix. One solves 
this as an eigenvalue equation of the form 



dqK(p,q;5)f(q) = A(S)f(p), (6.16) 

for a range of test values of 5 until an eigenvalue A (5) = 1 is obtained. 

The numerical solution to Eq. fl6.16| ) and its analogues for the negative parity states 
considered in Sections V and VI was carried out by linearly interpolating the functions 
contained in / on a 10 x 10 grid and then linearly interpolating the integrand on a 50 x 50 
grid. A cutoff of |q| , q^ < 10 in the dimensionless units defined by Eq. ( |6.9|) was found to 
be adequate. 

The calculations were done for the dimensionless parameter values jl = 3.5, 5.0 and 6.0. 
These values are chosen to give a broad range of jl while ensuring that the light propagator 
is confining (i.e. free from poles on the real timelike axis). Poles on the imaginary /c 4 -axis 
rendering the heavy the heavy quark propagator non-confining were sufficiently far away 
from the region of integration of the BSE to be considered irrelevant. However, one must 
also consider poles elsewhere in the the complex plane of both the light and heavy quark 
propagators, namely those listed in Table |, and the consequent restriction Eq. fl6.12|) on 
numerically accessible values of 5. To allow for as large a range of accessible 5 as possible, 
the momentum partitioning r] opt from Table | was used. 



In Fig. 3 we plot A(5) for the degenerate 0~/l~ mesons described by Eqs. (|Q ) and 
(or equivalently ( |5.24j ) and ( p.25| )), and the orbitally excited 1~ meson described by 
Eqs. ( |5.26| ) and ([5.271) . Also indicated in Fig. 3 are values of <5 max for each value of jl. As 
S approaches £ max it is clear that the numerical calculation breaks down as the integration 
routine is unable to cope with the singularity in the integrand. It is also clear that any 
reasonable extrapolation of the smooth part of the curves would not give a solution to 
A (5) = 1 for any of the curves in the range 6 < <5 max - It also appears that carrying through the 
calculations for higher values of jl would be unlikely to improve the situation. We therefore 
conclude that, if the model is taken in its simplest form, namely the combination of bare 
quark-gluon vertex and Gaussian gluon propagator in a Feynman-like gauge, spurious poles 
encountered in solutions to the DSE prevent the existence of solutions to the meson BSE. 
If the model is to be applied to the physics of D and B mesons, less crude approximations 
must be employed. 



VII. CONCLUSIONS. 



We have developed a non-perturbative treatment of heavy quarks which leads in a natural 
way to a dynamically generated self energy contribution to the heavy quark propagator. 
The treatment is based on a description of hadronic phenomena in terms of approximate 
Dyson-Schwinger equations which has proved highly successful within the light meson sector, 
namely the DSE technique j4|). 

The extension of the DSE technique to heavy quarks is non-trivial. The heavy quark 
self energy must be calculated near the bare current quark mass pole p 2 = -rnq, where p 
is the Euclidean momentum. This is accomplished by a change of momentum integration 
variable in the Dyson-Schwinger equation and an expansion in I/uiq. The dressed heavy 
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quark propagator is then used in conjunction with light quark propagators obtained from 
the conventional DSE technique to obtain a ladder approximation Bethe-Salpeter equation 
for D and B mesons to leading order in l/mq. The BSE gives the meson masses in the form 
rriQ + 5, where tuq is essentially infinite. At this level of approximation, the predictive power 
of the model lies not in calculating absolute meson masses, but in calculating mass differences 
between different spin states or different light quark content. We recover the degeneracy 
between pseudoscalar and vector mesons predicted by heavy quark effective theory M in the 
limit rriQ — > oo. 

We have carried out a numerical analysis of the model in the limit of infinite mass heavy 
quark, uiq — > oo, and massless light quark, m q = 0. Our calculation is not intended as 
a serious attempt to fit the D and B meson spectrum. Instead it is a feasibility study 
to determine to what extent the method is computationally viable. Solution of the BSE 
is numerically demanding, and the choice of model gluon propagator employed herein is 
primarily designed to reduce the numerical integration from three to two dimensions. We 
therefore use as input a particularly simple model gluon propagator, namely a gaussian 
propagator in a Feynman-like gauge. Because we are working in the chiral limit of the light 
quark, all dimensionful parameters can be scaled by the the width of the gaussian function, 
and the model is defined by a single adjustable parameter. 

We find the analytic structure of the heavy quark meson is consistent with that required 
for the change of integration variable used in the DSE. However we cannot unambiguously 
say whether the DSE leads to a confined heavy quark. Our bare vertex approximation to 
the DSE yields a succession of masslike poles which move off the timelike axis one by one as 
the parameter in the model gluon propagator is varied. As we point out below, we believe 
that these poles are a spurious artifact of the approximations employed. 

The combination of bare vertex and gaussian gluon propagator leads to poles in both the 
light and heavy quark propagators in the timelike half of the complex momentum plane over 
the broad range of input parameters studied. While poles off the real timelike momentum 
axis do not in principle spoil quark confinement, they can prevent successful solution of the 
BSE, which samples the quark propagators throughout a region of the complex plane. For 
the simple model set out above this is indeed the case. We have demonstrated numerically 
that the poles generated by the bare vertex approximation restrict the meson mass mq + 5 
to values of 5 which do not include any bound state solutions. 

In order to proceed further it will be necessary to develop an improved model which is 
not plagued by spurious singularities. In ref. |1| it is demonstrated that if the bare quark 
gluon vertex is replaced by a vertex ansatz satisfying the Ward-Takashi identity, such as the 
Ball-Chiu vertex P0| , the quark propagator becomes an entire function. This result holds 
independently of the current quark mass, and so must also hold in the heavy quark limit. 

Ideally one would like to incorporate quark propagator which is an entire function into 
our heavy quark formalism in a way consistent with successful DSE technique studies of light 
quark mesons. Of particular importance to the light meson sector is the requirement that the 
model should produce a Goldstone pion in the chiral limit. This is achieved by ensuring that 
the choice of Bethe-Salpeter kernel is compatible with the choice of quark-gluon vertex in the 
DSE [|18j. While there is no practical way at present to achieve such a compatibility with the 
Ball-Chiu vertex, good results can be achieved for the light pseudoscalar and vector meson 
octets by inverting a bare vertex approximation DSE from an entire light quark propagator 
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to construct a separable approximation to the Bethe-Salpeter kernel Extension of the 

separable approximation in a similar manner to the heavy quark sector, beginning with the 
light quark propagators employed in ref. || and a heavy quark propagator obtained from 
the Ball-Chiu vertex is perhaps the most promising way to proceed at this point. One 
potentially then has a single non-perturbative model of meson dynamics applicable to the 
entire range of meson masses. 

Finally we note that, if the DSE technique is taken as a serious description of hadronic 
phenomena, we are led unavoidably to a treatment generically of the kind set out in this pa- 
per. The above analysis illustrates for the first time the importance of the analytic structure 
of heavy quark propagators and its significance in determining the heavy meson spectrum. 
Our analysis leads to a heavy quark self energy which will contribute as significantly to 
the dynamics of the D and B mesons as that of the light quarks, and therefore cannot be 
ignored. Furthermore we are unaware of any existing treatment of the heavy quark self 
energy which is genuinely non-perturbative and has the potential to provide a qualitative 
explanation of heavy quark confinement. For these reasons we believe it is important to 
persevere with the application of the DSE technique to the heavy quark sector. 
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APPENDIX. 



We derive here a set of coupled integral equations from the DSE for the heavy quark 
propagator to order 1/rriQ. From Eqs. Q2.1D, ( |2.3| ) and ( |2.4j ) we obtain 



S(p) 



1+ 



1 



+ 



1 



2m. 



Q 



2 ik ■ v + Y, B — Y>a 



ik ■ v + Tj B — Y, A 2 (ik ■ v + T, B — a^a) 
where, on choosing = (1, 0), 

£4 = X4O4, |k|), E B = E B (k^ |k|). 
Clearly and E# should depend only on one variable, namely 



+ O (m' c 



V 



-m 2 + 2imnk 4 + /c 4 2 + |k| . 



(Al) 



(A.2) 



It follows that 
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E A (fc 4 , |k|) = £ A ^ 4 - ^— |k| 2 + O (m Q 2 ) , 0^ 

% ,,,2^(^4,0) 



2mr 



<9fc 4 



+ o K 2 ) • 



We now make the expansion 



1 



Then 



E^O) = E° (fc 4 ) + E^ 4 ) + 

m <2 



£ A (fc 4 , |k|) = £° (fc 4 ) + — [Ei(*0 - - |k| 2 E°'J + O (m Q 2 
= S o (fc4 ) + _L^(A;4,|k|) + 0(m Q 2 ), 



(A.3) 



(A.4) 



(A.5) 



and similarly for Eg. (E^ is only defined for convenience here. Eventually we shall be 
solving for E° , T,° B , T,\ and E#, which only depend on fc 4 .) For the propag ator Eq. flPD 
we then have 



Sip) 



+ 



1 



2 iA; 4 + - E^ 2m 



1+ ^ 4 £° + (E^) 2 - (Ej ) 2 + fc 4 2 + \k\' + 2 (g, - Ej 

(ifc 4 + - EQ )2 

The left hand side of the DSE Eq. (|2.2j ) then becomes 



+ O (mj) 



(A.6) 



m Q 



■E. + E, 



1 



^° A (h) + Z B (k A ) + % ^° A (k A )- ^\(k A , |k|) + t B (h, |k|) 

m Q 1 



(A.7) 



On the right hand side, we substitute in Eq. (|A.6|) and change the variable of integration to 
k' defined by g M = imqv^ + k' . Then setting k^ = (/c 4 , 0), and projecting out coefficients of 
tf, ft and I gives the following set of integral equations: 



4 r d A k' 



(27T) 



ik 4 ~\~ E 1 ^/^) E^(/c4) 



(A.8) 



(A.9) 



4^ 4 £° +3(£ A ) 2 + (fc 4 ) 2 + |k'| 



ik'd -\- E^ 



2 ( Tj B — T, A 



IlkfE . 



(ik> 4 + £ A )i 



(A. 10) 
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2E« 



ik^ + E^ 



4i^E° + 3(E° ) 2 + (^) 2 + |kf + 2 (Ek - Y,\ - \ |kf E° ') 



(A.11) 



Once a model quark propagator A(fc — fc') is specified, these equations can, in principle, be 
solved numerically to give the heavy quark propagator to first order in 1/toq. 
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TABLES 



TABLE I. Position of the poles sP ole in the light quark propagator and k^ ole in the heavy quark 
propagator, and the maximum accessible values of <5 max and corresponding momentum partitioning 
parameter r] opt - Dimensionful quantities can be extracted from the dimensionless quantities listed 
via Eqs. Q to Q. 





-pole 


£,Poic 
"'4 


^max 


Vopt 


2.5 


-0.126 + 0.000* 


1.802 + 0.957* 


1.307 


0.268 


3.0 


-0.540 + 0.000* 


2.285 + 1.034* 


1.769 


0.415 


3.5 


-1.829 + 0.919*' 


2.769 + 1.102*' 


2.494 


0.558 


4.0 


-1.646 + 2.234*' 


3.255 + 1.162*' 


2.649 


0.561 


5.0 


-1.246 + 4.546*' 


4.229 + 1.264*' 


2.991 


0.577 


6.0 


-0.327 + 7.216* 


5.207 + 1.351*' 


3.294 


0.589 


7.5 


2.132 + 11.711« 


4.680 + 2.649*' 


4.859 


0.455 
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FIGURE CAPTIONS. 

Figure 1: The heavy quark self energy obtained by solving Eq. (|6.3|): ReS(fc 4 ) (dashed 
curve) and — ImS^) (solid curve) for fi = n/^/a = 3.5. 

Figure 2: The position k\ ole of the poles in the heavy quark propagator for the values of fi 
indicated. Only poles for which Re k 4 > are shown, though the poles come in pairs 

j pole j f pole * 

k\ and — k% 

Figure 3: Eigenvalues A(<5) from Eq. fl6.16|) for the degenerate 0~/l~ states (solid curves) 
and the orbitally excited 1~ state (dashed curve) calculated with gluon propagator 
parameter fi = 3.5 (triangles), 5.5 (circles) and 6.0 (squares). Also shown as vertical 
dotted lines are the numerical limits <5 max for each value of fi. 
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